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WHAT IS CLAIMED IS: 

1. A method for correcting the color contrast of a scene, the scene including an 
intensity spectrum at each of a plurality of pixels, the method comprising the steps 
of: 

a) providing a red image, a green image, and a blue image, each image having a 
pixel value at each of the plurality of pixels; 

b) computing a center red response, a center green response and a center blue 
response based on said images; 

c) computing a surround red response, a surround green response and a surround 
yellow response based on said images; 

d) computing a red, a green and a blue on-center opponent and filtered opponent 
response, based on said center and surround responses; 

e) computing a red, a green and a yellow off-center opponent and filtered 
opponent response based on said center and surround responses; 

f) computing a red, a green and a blue double-opponent response and a 
corresponding filtered double-opponent response based on said on-center and 
off-center filtered opponent responses; 

g) computing a red, a green and a blue do-remote signal based on a set of 
responses selected from the group consisting of said on-center filtered opponent 
responses and said filtered double-opponent responses; and 

h) for each pixel: correcting each of said red, green, and blue double-opponent 
responses for color contrast using respectively said red, green and blue do- 
remote signals, thereby producing corrected red, green and blue double- 
opponent responses. 
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2. The method of claim 1, further comprising the steps of: at each pixel: 

i) computing a yellow center response and a blue surround response; 

j) computing a yellow double-opponent response and a corresponding yellow 

filtered double-opponent response based on said yellow center and said 

blue surround responses; 
k) computing a yellow do-remote signal based on a response selected from 

the group consisting of said yellow center response and said filtered 

yellow double-opponent response; and 
1) for each pixel: correcting said yellow double-opponent response for 

color contrast using said yellow do-remote signal, thereby producing a 

corrected yellow double-opponent response, 

3. The method of claim 1, wherein said step of providing a red image, a green 
image, and a blue image includes: at each pixel: 

i) multiplying the intensity spectrum by a spectral response function of a 
red photoreceptor, thereby providing a red spectral product; 

ii) multiplying the intensity spectrum by a spectral response function of a 
green photoreceptor, thereby providing a green spectral product; 

iii) multiplying the intensity spectrum by a spectral response function of a 
blue photoreceptor, thereby providing a blue spectral product; 

iv) integrating said red spectral product, 

v) integrating said green spectral product, and 

vi) integrating said blue spectral product. 
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4. The method of claim 1, wherein said step of computing a center red response, a 
center green response and a center blue response includes, for each said image, 
convolving said image with a center local spatial filter; and wherein said step of 
computing a surround red response, a surround green response and a surround 
yellow response includes, for each said image, convolving said image with a 
surround local spatial filter. 

5. The method of claim 2, wherein said step of computing said red on-center 
filtered opponent response includes subtracting said green surround response 
from said red center response, said step of computing said green on-center 
filtered opponent response includes subtracting said red surround response from 
said green center response, said step of computing said blue on-center response 
includes subtracting said yellow surround response from said blue center 
response. 

6. The method of claim 2, wherein said step of computing said red off-center 
filtered opponent response includes subtracting said green surround response 
from said red center response, said step of computing said green off-center 
filtered opponent response includes subtracting said red surround response from 
said green center response, and said step of computing said yellow off-center 
filtered opponent response includes subtracting said yellow surround response 
from said blue center response. 
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The method of claim 1, wherein said step of computing said red, green and blue 
filtered double-opponent responses based on said on-center and off-center 
filtered opponent responses further includes computing respective red, green 
and blue double-opponent center responses and respective red, green and blue 
double-opponent surround responses. 

The method of claim 7, wherein said step of computing each said double- 
opponent center response includes convolving a center filtered response with a 
center spatial weight function 

The method of claim 1, wherein said step of computing each said double- 
opponent surround response includes convolving a surround filtered response 
with a surround spatial weight function 

The method of claim 8, wherein said center spatial weight function is an 
exponentially decaying function. 

The method of claim 9, wherein said surround spatial weight function is an 
exponentially decaying function. 

The method of claim 7, wherein said step of computing each said do-remote 
response includes convolving a corresponding response selected from the 
group consisting of a said on-center filtered opponent response and a said 
double-opponent center response, with a remote spatial weight function. 
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13. The method of claim 12, wherein said spatial weight function in chosen from 
the group consisting of exponentially decaying functions and Gaussian 
functions. 

14. The method of claim 1, wherein said step of correcting each of said red, green, 
and blue double-opponent responses for color contrast includes the steps of: 
for each said double-opponent response 

a) computing a respective adaptive function Gb ; and 

b) computing a respective adaptation factor, based on said respective 
adaptive function. 

15. The method of claim 14, further including applying said respective 
adaptation factor as a respective semi-saturation factor in a Naka-Rushton 
type equation operating on each of said double-opponent responses, 

16. The method of claim 2, wherein said step of correcting said yellow double- 

opponent response for color contrast includes the steps of: for each said 
double-opponent response 

c) computing an adaptive function Gb ; 

d) computing an adaptation factor, based on said adaptive function. 

17. The method of claim 14, further including applying said adaptation factor as 
a semi-saturation factor in a Naka-Rushton type equation operating on said 
yellow double-opponent response. 
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18. The method of claim 15, wherein said adaptive function is identical to said do- 
remote response. 

1 9. The method of claim 15, wherein said adaptive function is time-dependent. 

20. The method of claim 17, wherein said adaptive function is identical to said do- 
remote response. 

2 1 . The method of claim 1 7, wherein said adaptive function is time-dependent. 

22. The method of claim 19, wherein said time-dependency of said adaptive 
function is obtained by a convolution of said double-opponent response with a 
temporal filter. 

23. The method of claim 21, wherein said time-dependency of said adaptive 
function is obtained by a convolution of said double-opponent response with a 
temporal filter. 

24. The method of claim 22, wherein said temporal filter is an exponentially 
decaying time function. 

25. The method of claim 23, wherein said temporal filter is an exponentially 
decaying time function. 
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The method of claim 22, wherein said temporal filter is normalized by a 
denominator equal to x b (t) = x m /( 1 + abs( G(t) - G b (t) )/G„ ) where x m is an 
expected maximum value of said decaying time function, and where G n is a 
normalization constant. 

The method of claim 23, wherein said temporal filter is normalized by a 
denominator equal to x b (t) = x m /( 1 + abs( G(t) - G b (t) )/G n ) where x m is an 
expected maximum value of said decaying time function, and where G n is a 
normalization constant. 

The method of claim 1, further comprising the step of: inversely transforming 
said corrected double-opponent responses into corresponding new red, green 
and blue center responses. 

The method of claim 28, wherein said inversely transforming includes 
transforming said new double-opponent cell responses into new opponent cell 
responses. 

The method of claim 2, further comprising the step of: inversely transforming 
said corrected yellow double-opponent response to obtain a new yellow 
opponent response. 

A method for adjusting an achromatic contrast of a scene, the scene including 
an intensity spectrum at each of a plurality of pixels, the method comprising 
the steps of: 
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a) providing an image that has an intensity value at each of the plurality 

of pixels; 

b) obtaining an adapted opponent center response using a plurality of said 
pixel intensity values, and 

c) at each pixel, correcting the achromatic contrast using said adapted 
opponent center response. 

32. The method of claim 3 1 , wherein said step of obtaining an adapted opponent 
center response includes: 

i) calculating an opponent center response; 

ii) providing a center adaptation factor, and 

iii) combining said opponent center response and said center adaptation 
factor. 

33. The method of claim 32, wherein said adaptation factor includes a remote 
center adaptation term. 

34. The method of claim 33, wherein said adaptation factor further includes a 
local center adaptation term. 

35. The method of claim 32, further comprising obtaining an adapted opponent 
surround response, wherein said step of correcting for achromatic intensity 
contrast includes subtracting said adapted opponent surround response from 
said adapted opponent center response. 
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36. The method of claim 35, wherein said step of obtaining an adapted opponent 
surround response includes: 

i) calculating an opponent surround response; 

ii) providing a surround adaptation factor, and 

iii) combining said opponent surround response and said surround 
adaptation factor. 

37. The method of claim 36, wherein said surround adaptation factor includes a 
remote surround adaptation term. 

38. The method of claim 36, wherein said surround adaptation factor further 
includes a local surround adaptation term. 

39. The method of claim 3 1 , wherein said step of providing pixel intensity values 
includes: at each pixel : 

i) multiplying the intensity spectrum by a spectral response function, 

thereby providing a spectral product; and 

ii) integrating said spectral product. 

40. The method of claim 31, wherein said step of obtaining an opponent center 
response includes convolving each said pixel intensity value with a center 
spatial weight function. 
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41. The method of claim 36, wherein said step of obtaining an opponent surround 
response includes convolving each said pixel intensity value with a surround 
spatial weight function. 



42. A method for adjusting an achromatic contrast of a scene, the scene including 
an intensity spectrum at each of a plurality of pixels, the method comprising 
the steps of: 

a) providing an image having a pixel value at each of the plurality of pixels; 

b) computing a center response based on said image; 

c) computing a surround response based on said image; 

d) computing an on-center filtered opponent response, based on said center and 
surround responses; 

e) computing an off-center filtered opponent response based on said center and 
surround responses; 

f) computing a double-opponent response and a corresponding filtered double- 
opponent response based on said on-center and off-center filtered opponent 
responses; 

g) computing a do-remote response; and 

h) for each pixel: correcting said double-opponent responses for achromatic 
contrast using said do-remote signal. 



43. The method of claim 42, wherein said step of providing an image includes: at 
each pixel: 

i) multiplying the intensity spectrum by a spectral response 
function to obtain a spectral product; and 
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ii) integrating said spectral product. 

The method of claim 42, wherein said step of computing a center response 
includes convolving said image with a center local spatial filter. 

The method of claim 42, wherein said step of computing said on-center filtered 
opponent response includes subtracting an on-center surround response from an 
on-center center response. 

The method of claim 42, wherein said step of computing said off-center filtered 
opponent response includes subtracting an off-center surround response from an 
off-center center response. 

The method of claim 42, wherein said step of computing said double-opponent 
response and said filtered double-opponent response includes computing a 
double-opponent center response. 

The method of claim 42, wherein said step of computing said double-opponent 
response and said filtered double-opponent response further includes computing 
a double-opponent surround response. 
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